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Abstract 

Geometrical random multiplicative cascade processes are often used to model positive- 
valued multifractal fields such as the energy dissipation in fully developed turbu- 
lence. We propose a dynamical generalization describing the energy dissipation in 
terms of a continuous and homogeneous stochastic field in one space and one time 
dimension. In the model, correlations originate in the overlap of the respective space- 
time histories of field amplitudes. The theoretical two- and three-point correlation 
functions are found to be in good agreement with their equal-time counterparts 
extracted from wind tunnel turbulent shear flow data. 
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Whenever strongly anomalous, intermittent fluctuations, long-range correla- 
tions, multi-scale structuring and selfsimilarity go hand in hand, the label 
'multifractality' is attached to the underlying process. While in this Letter 
we have in mind fully developed turbulence of fluid mechanics [1], such pro- 
cesses occur in various other fields such as formation of cloud and rain fields 
in geophysics [2], internet traffic of communication network engineering [3], 
and stock returns in finance. [4], to name but a few. 

Random multiplicative cascade models (RMCMs) are commonly used to model 
and visualize such phenomena since they generally exhibit multifractality and 
reproduce the abovementioned properties [5]. They are usually implemented 
through a scale-independent cascade generator which produces a nested hier- 
archy of scales and multiplicatively redistributes the local measure. 

In fully developed turbulence, RMCMs have often been employed to model 
the energy flux through inertial-range scales. Due to their multiplicative na- 
ture, they can easily reproduce multifractal scaling exponents associated with 
the energy dissipation [6] , the latter representing the intermittency corrections 
[1]. Although the link between such models and the Navier-Stokes equation re- 
mains unclear, recent investigations on multiplier distributions [7,8] and scale 
correlations [9] have shown that RMCMs do appear to contain more truth 
than might reasonably be expected from their phenomenological basis. 

Nevertheless, these discrete RMCMs are purely geometrical constructs and in- 
capable of describing causal dynamical effects of the turbulent energy cascade. 
A generalization in this direction is clearly desirable. Hence, in this Letter, we 
present a dynamical RMCM in 1+1 space-time dimensions which respects 
causality and homogeneity. It is related to recent, related efforts [10,11,12,13], 
but goes beyond them in its dynamical interpretation. We first show how this 
model yields multifractal scaling for arbitrary n-point correlation functions, 
proceeding thereafter to compare equal-time two- and three-point correlation 
functions to their counterparts obtained from wind-tunnel turbulent shear flow 
data. 

Our dynamical RMCM is constructed by analogy to the geometrical case, in 
which the amplitude of the positive-valued energy-dissipation field, resolved 
at the dissipation scale i], is defined as the product of independently and 
identically distributed random weights q(lj), 



where lj is an element of a nested hierarchy of scales rj — lj < lj — L/A- 7 < 
l — L with 0<j<J the "cascade generation" and A>1 the discrete scale step. 
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The integral length L and the dissipation length rj represent, respectively, 
the largest and smallest length scale of the process. The geometrical RMCM 
furthermore requires (q) = 1 because of conservation of energy flux. 

We generalize (1) by assuming that e is again the multiplicative product of a 
stochastic field, but that this field is now defined on continuous 1+1 spacetime: 

(oo oo *\ 

J dt' J dx'f(x -x',t- t')7<y , | , (2) 



where / is the "index function" described below and by assumption 7(2;, t) ~ 
S a ((dxdt) a _1 <r, — is a Levy-stable white-noise field with index 0<a<2 
[14]. For a— 2, this corresponds to a non-centered Gaussian white-noise field. 
From its characteristic function, (expjrry}) = exp{ — (a a n a )/(cos(7ro;/2)) + 
fin} with a^l, the parameter \i is fixed to /i — cr a / cos(7ro;/2) in order to 
satisfy the requirement (exp{7}) = 1. 

Causality, i.e. the requirement that e(x,t) depends on the past but not the 
future, dictates that the index function f(x — x',t — t') must be zero for 
t — t'<0. Demanding also spatial symmetry around x, we are led to the form 



f(x-x',t-t') = < 



1 ( 0<t - f<T, -g(t - t')<x - x'<g(t -t')), 
(otherwise) . 



As illustrated in Fig. 1, the causality cone g(t — t') incorporates a correlation 
time T and a correlation length L with g(T) = L/2. The exponent of the 
ansatz (2) can be thought of as a moving average over the stable white-noise 
field. 

According to (3), the time integration in (2) runs over 0<t — t'<T. Since to 
any given time tj there corresponds a length scale lj—2g(tj), there is a joint 
hierarchy of length and time scales, so that (2) factorizes into integrals of 7 
over the separate slices shown in Fig. 1, 



t-u 



q(lj) = exp < J dt' J dx'"f(x',t' 

t-tj-i x-g(t-t') 



(4) 



In order to interpret q(lj) as a random multiplicative weight, its probability 
density needs to be independent of scale. Since the , y(x',t') are i.i.d., the 
integration domain of (4) must therefore be independent of the scale index j. 
Together with the the boundary conditions g(T—AT L )=L/2 and g(AT v )=r]/2, 
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this fixes the causality cone to 



g[l ' ~ l \ ( L -V) (T-AT L -(t-t')) y°) 



r, (T-AT L -AT n ) 



for times AT v <t — t'<T — AT^. To complete the picture, we need to specify 
g(t - t') for < t - t' < AT V and T - AT L <t-t'<T. Since on physical 
grounds we expect AT^T, the simplest choice is AT V = 0. For the remaining 
parameter AT L) we assume AT L <C T; it will be specified more fully below. 

The construction proposed in Eqs. (2)-(5) guarantees that the one-point statis- 
tics of the dynamical RMCM are identical to its geometrical counterpart. In 
order to qualify for a complete dynamical generalization, not only the one- 
point statistics but the n-point statistics in general should match. Hence, we 
now consider the equal-time two-point correlator with t\=t 2 =t and Ax = 
x 2 -x l > 0, 



P , A ,_ (e^(x 1 ,t)e^(x 2 ,t)) 



(e n ^X!,t)) (e n *(x 2 ,t)) 

{D(Ax) ni+n ' 2 } 
(D(Ax) n i) (D(Ax) n i) ' 



(6) 



As shown in Fig. 2a, the correlation between two points a distance 77 < Ax < L 
apart stems from the overlap region of the two index functions f(x' — x\,t) 
and f(x 2 — x', t). This explains why, in the second step of (6), R can be written 
solely in terms of integrals D over the overlap region, 

/t- 9 ^){Ax/2) x 1+ g(t-t') \ 

D(Ax) = expl J dt' J dx'-f(x',t') , (7) 

\ t-T x 2 -g(t-t') J 



since the contributions from the non-overlapping regions are statistically in- 
dependent and hence factorize and cancel. Introducing the spatio-temporal 
overlap volume 



t-g(~V (Ax/2) X1 +g(t-t') 

V(Ax)= J dt' J dx' 

t-T X2 -g( t -t') 

_ r,L{T - AT,) / L \ fv(T - AT VJ 



- (L- n) to lsJ-i^#- AIi J (L - A * ) ' (8) 

and employing basic properties of stable distributions [14], the expectation of 
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the n-th power of D is found to be 

(D(Axr) = exp [^aV(Ax)(n - . (9) 

Defining the multifractal scaling exponents r{n) = r(2)(n— n Q )/(2— 2 Q ), with 
r(2) = (a a /cos^f)(2-2 a )r]L(T-AT v )/(L-r]), as well as T[n u n 2 ] = r(ni + 
n 2 ) — ririi) — r(n 2 ), substitution of (9) into (6) leads to the final expression 
for the equal-time two-point correlator: 



Rni,n2 {Ax) 



exp 



-r[n 1 ,n 2 \ 1 




(10) 



Equal-time two-point statistics of our dynamical RMCM in 1+1 dimensions 
hence show multiscaling behavior for r/<Ax<^L, in complete analogy to the 
findings of the corresponding geometrical RMCM [15]. We also note that set- 
ting ATl = (T — AT v )r)/(L — r/) eliminates the second factor in (10) leaving 
R ni , n2 (Ax) = (L/ Ai) 7 '" 1 '" 2 ' to scale rigorously. 

Turning to temporal two-point correlations, we follow the same recipe as for 
the above. Correlations in this case arise from the overlap volume illustrated 
in Fig. 2b, and lead, after an analogous straightforward calculation, to 



(e n ^x,h)) (e n *(x,t 2 )) 

At - AT L n / At- ^ TL ^-r[ni,n 2 ](l-AT L /T) 



T - AT L L V T — AT L 

with At = t 2 — t\ and xi=x 2 =x. For simplicity, the parameter AT V has been 
set to zero. The last step of (11), valid for AT L <tiAt<T and r\^L only, shows 
that the temporal two-point correlator has scaling exponents identical to those 
of the equal-time counterpart. 

Although respective spatio-temporal overlap volumes are more complicated, 
two-point spacetime correlations with both Ax ^ and At ^ can also be 
derived; see Ref. [16] for a complete analysis. Here, we prefer to continue with 
equal-time three-point correlations; their generalization to equal-time n-point 
correlations is straightforward and explicit expressions can again be found in 
Ref. [16]. The corresponding overlap volumes, illustrated in Fig. 2c, represent 
the starting point for a calculation analogous to (7)-(9), which leads to 
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Rn\,n2,ns \£\ j %2i -£3) 

(s ni (x 1 ,t)e n2 (x 2 ,t)e n3 (x 3 ,t)) 
~ (e n i( Xl ,t)) (e n i(x 2 ,t)) (e n *(x 3 ,t)) 

\ r[ni+n2,n3]-r[n2,n3] , ^ 



X 3 - Xi 

exp < - 



x 2 - X 1 



r[ra 1 ,n 2 J 



L 



x 3 - x 2 



r[n 2 ,nz\ 



+r[ni,n 2 ] 1 



(L-t,)ATl 
V(T-AT V ) / 
x 2 -x l 



r T[ni+n 2 ,n 3 }-T[n 2 ,n 3 ]) ( 1- 



+ r[n 2 ,n 3 ] 1 



0]} 



(12) 



with x\ < x 2 < x 3 and r\ < x^ — Xj < L for all i > j = 1,2,3. In the case of 
n\—n 2 — n 3 — \ and small separations \xi — Xj\ <C L, or for all separations 
if the parameter ATl is fine-tuned, this simplifies to 



-Ri i i,i(x 2 — x 1 =const, 0:3=0;) 

v r(2) / T xr(2) 



(-^-T' P-) n ' {x l<X < x 2 ) 



The comparative ease with which the three-point expressions (12) and (13) 
were derived can be traced back to the fact that the present model incor- 
porates spatio-temporal homogeneity from the very beginning. This is to be 
contrasted with the geometrical RMCMs of (1) which, due to their hierarchical 
structure, are not translationally invariant in space. This non-invariance feeds 
through to all n-point observables and has to be removed at considerable cost 
through successive spatial sampling [17] before the latter can be compared to 
experimental counterparts. 

While we expect the model proposed in this Letter to find application in many 
different phenomena, we demonstrate its qualities through comparison with 
fully developed turbulence data. A velocity record of the longitudinal compo- 
nent, obtained in a wind-tunnel shear flow experiment [18], was transformed 
into a one-dimensional spatial record of the positive- valued surrogate energy 
dissipation e(x) = I5v(d x v x ) 2 , where v is the viscosity. The sampled two-point 
correlations of order ri\ — n 2 — 1 and ri\ = 2, n 2 = 1 are plotted in Fig. 3. 
Within the inertial range 77 \x 2 — xi\ C L, the data reveals rigorous power- 
law scaling with exponents r[l, 1] = 0.184 and r[2, 1] = 0.34. This fixes the 
intermittency exponent r(2) = 0.184 and the stable index a = 1.80, which 
are the relevant model parameters for multifractal scaling. Once these have 
been fixed, no further room for adjustment is left for the theoretical three- 
point correlation (13), which is compared to its experimental counterpart in 
Fig. 4. Independent of the various combinations for the two-point distances 
r] <^ \xi — Xj\ <C L, the agreement between model and data is remarkable. 
This demonstrates that the proposed stochastic process, whose parameters 
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have been fixed from lowest-order two-point correlations, is capable to describe 
the equal-time multivariate statistics of the turbulent energy dissipation be- 
yond two-point order. It proofs also that the turbulent energy cascade can be 
thought of as a consistent multifractal process. 

The dynamical RMCM presented here is a generalization of the geometrical 
RMCM. By construction, it is causal, continuous and homogeneous, does not 
make use of a discrete hierarchy of scales and stochastically evolves a positive- 
valued field in one space and one time dimension. Several generalizations of this 
new model come to mind immediately, such as the stochastic evolution in n+1 
dimensions with the optional inclusion of spatial anisotropy, the use of inde- 
pendently scattered random measures to describe deviation from log-stability 
[19], the discretisation of space-time into smallest cells to model dissipation, 
and a dynamical RMCM for vector fields to model the turbulent velocity field. 
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Fig. 1. Causal space-time "cone" for the positive-valued multifractal field e(x,t). 
All field amplitudes j(x',t') inside the causal space-time "cone" bordered by the 
index function (3) contribute multiplicatively to e. See text for further detail. 




(«) (&) to 



Fig. 2. Spatio-temporal overlap volumes (shaded) producing the correlation for 
the (a) equal-time and (b) temporal two-point correlator, as well as (c) for the 
equal-time three-point correlator. To simplify visualization, parameters have been 
set to AT V = AT L = 0. 
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Fig. 3. Two-point correlator (e ni (x 1 ,t)e n2 (x 2 ,t)) /((e ni (x 1 ,t)}(e n2 (x 2 ,t))) of or- 
ders (a) ni=ri2=l and (b) ni=2, n 2 =l for the experimental shear-flow dataset [18] 
with Taylor Reynolds number R\ = 860 and integral length scale L = 19607? as a 
function of the distance \x 2 — x±\ in units of the dissipation length rj. The insets 
represent the compensated plots, where the two-point correlators have been divided 
by (\x 2 - xi|/r/)- T t ni ' n2 ] with r[l, 1] = 0.184 and r[2, 1] = 0.34, respectively. 
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Fig. 4. Comparison between expression (13) and the experimentally extracted 
three-point correlator {e{x\, t)e(x2, t)e(xs, t))/{e(x, t)) 3 . The same data set as in 
Fig. 3 has been used. The scaling exponents r[l, 1] = 0.184 and r[2, 1] = 0.34 
have already been fixed by the two-point correlators. In part (a), the experimental 
three-point correlator is shown for a fixed two-point distance X2~ x\ = 336ry. The two 
insets represent the comparison with (13) for the two regimes rj < x^—xi < X2—X1 
and X3— x\ > X2~x\. Part (b) focuses on the regime 77 < X2— x\ < X3—X1 with 
fixed (23— xi)/rj = 336, 672 and 1344, respectively, while (c) focuses on the regime 
^3—^1 > X2—X1 with fixed (x2—xi)/rj = 84, 168 and 336, respectively. For clarity, 
from left to right the curves of (c) have been shifted by 0, 0.4 and 0.8 along the 
y-axis. 
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